Efficient and robust initialization of a qubit register with fermionic atoms 
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We show that fermionic atoms have crucial advantages over bosonic atoms in terms of loading in 
optical lattices for use as a possible quantum computation device. After analyzing the change in 
the level structure of a non-uniform confining potential as a periodic potential is superimposed to 
it, we show how this structure combined with the Pauli principle and fermion degeneracy can be 
exploited to create unit occupancy of the lattice sites with very high efficiency. 



Ultra-cold neutral atoms list among the most promis- 
ing candidates for quantum computing (QC) applica- 
tions 0. Several possible QC schemes using cold atoms 
have recently been proposed (H, and various groups are 
presently working for their experimental realization 
The common starting point is the creation of a quantum 
register, i.e. a scalable array of well-characterized qubits 
The qubits are usually indentified with two internal 
states, typically hyperfine states, of a single atom. The 
atom-qubits need to be spatially well separated, and yet 
to be at a distance which allows to efficiently activate 
their interaction. These requirements are nearly ideally 
fulfilled by a periodic potential with exactly one atom lo- 
calized at each lattice site. Such a periodic potential can 
be generated, for instance, by counter-propagating laser 
beams or with microscopic magnetic surface traps (atom 
chips). As a natural consequence of the vast knowl- 
edge developed in the last decade in the field of Bose- 
Einstein condensation Q , until now most efforts towards 
atomic QC implementations have been directed at the 
use of bosonic atoms. Only recently an increasing inter- 
est has been devoted to the cooling and manipulation of 
fermionic atoms in confining potentials Q and opti- 
cal lattices 8]. In particular the possibility to create a 
quantum register with fermionic atoms in a homogeneous 
optical lattice via coherent filtering has been proposed in 

i 

In a uniform system, a lattice with a single atom per 
site can be created by confining N atoms in a one- 
dimensional box of length L, and by subsequently su- 
perimposing a periodic potential of spatial period d, so 
that the number of lattice sites in the box is L/d = N. If 
the periodic potential is raised sufficiently slowly and if 
the effective atom-atom interaction is repulsive, the sys- 
tem ends up in its ground state which has an average 
occupation of one atom per site. One could equivalently 
use bosons with a positive scattering length or polarized 
fermions, in which case the Pauli principle provides a 
strong effective repulsion. In either case, however, this 
simple scheme has significant problems. To avoid the cre- 
ation of defects like holes or multiply occupied wells, it 
requires equal number of particles and lattice wells be- 
sides a long adiabaticity time. 

In this Letter we show that, in the presence of a non- 
uniform external confinement instead, fermionic atoms 



have crucial advantages over bosonic atoms and we pro- 
pose a simple protocol to create a quantum register which 
directly exploits the pecidiar properties of the Pauli prin- 
ciple. 

Before discussing our proposal, it is useful to further 
consider the case of non-interacting fermions confined by 
a homogeneous box. In the absence of the periodic poten- 
tial, the A'^ fermions occupy the N lowest single-particle 
levels up to the Fermi energy ep = 7r^?i^A^/2mL^, with 
TO being the atomic mass. When the periodic potential is 
raised, in order for the number of lattice sites to be equal 
to the number of particles we should have = ep, where 
Er — TT^h^ /2md^ is the recoil energy of the lattice. In- 
creasing the depth of the optical lattice slowly enough to 
avoid excitations, the lowest A'' levels gradually cluster 
to form the first Bloch band, yet conserving their initial 
unit occupancy due to the Pauli principle. Eventually, in 
the limit s ^ oo the A^ levels become degenerate and the 
Wannier wavefunctions, which are N linear combinations 
of the Bloch states localized one in each well, are exact 
eigenstates of the system. Under realistic conditions this 
procedure has the following drawbacks: 

1. For finite s, the Wannier states are only approximate 
eigenstates and delocalize on a time scale of the order of 
the inverse of the bandwidth 6. This limits the time over 
which one can perform QC operations. 

2. Randomly distributed defects appear in the register 

(i) if ep ^ Er — the condition ep < Er results in empty 
sites, while ep > Er results in multiply occupied sites; 

(ii) if the initial temperature is not much smaller than 
the energy level spacing; 

(iii) if the adiabaticity condition is not perfectly fulfilled, 
i.e. if the raising time of the lattice is not large as com- 
pared to the inverse level spacing. 

The conditions required to create a lattice with unit 
filling along the lines specified above are unrealistically 
strict, and more complex schemes (see e.g. a^re 
needed in order to eliminate the defects produced in the 
register. 

When the Pauli principle is combined with a non- 
uniform confinement, the scenario changes completely. 
We first consider a one-dimensional geometry, which can 
be realized in practice with a tight radial confinement, so 
as to freeze excitations of the radial degrees of freedom. 
The extension to higher dimensions will be briefly dis- 
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cussed later. In the presence of an external confinement, 
Vext, the total trapping potential is 



5=7.2 



5=72 



V{z) = Vea:tiz) + sErSVC?{TTz/d), 



(1) 



where sEr is the lattice height. In order to present the 
corresponding loading procedure, we need to analyze the 
behavior of the energy eigenvalues and the corresponding 
eigenstates as a function of the optical potential depth. 
In Fig.^ we consider the specific case of a harmonic ex- 
ternal confinement Vext{z) — muj1[z — zo)^/2, where zo ( 
is an arbitrary offset between the center of the confining 
potential and the optical lattice. As the periodic poten- ^^^j 
tial is raised, the eigenstates change from the well known 
oscillator states at s = 0, to states localized at the lattice 
sites for large s. 
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FIG. 1: Energy levels relative to the ground state £"0 as a 
function of the optical potential depth s, for a harmonic con- 
finement with hoj = i5r/(27r^) and zo = d/3. 



It is a general feature that the eigenvalues of the initial 
external confinement (s = 0) can be grouped into two 
sets: 

(I) The states with energy smaller than the recoil energy 
change smoothly as s increases; 

(II) The levels with energy larger than the recoil energy 
show a more complicated behavior characterized by sev- 
eral avoided crossings. 



The reason for this distinction, and the mentioned lo- 
calization of the eigenstates, can be understood by look- 
ing at the behaviour of the eigenvalues and the eigen- 
functions, illustrated in Fig. [3 The figure shows the 
spatial position and the extension of the eigenfunctions 
versus the corresponding eigenvalues for two values of s. 
When s increases the eigenfunctions become more and 



FIG. 2: Spatial position and extension of the eigenfunctions 
versus the corresponding eigenvalues for s = 1.2 (a), and s — 
12 (b); Density profiles calculated for a number of atoms such 
that ep = Er at s ^ 0, for s = 1.2 (c), and s = 12 (d). The 
confining potential is the same as in Fig. Q 

more localized in space and a band-like structures ap- 
pear. The gaps between pseudobands open at the center 
of the trap, and the one between the first and the sec- 
ond band separates the initial eigenstates with energy 
smaller and larger than Er, i.e. the states in sets (I) and 
(II) above. The eigenfunctions with initial energy smaller 
than Er gradually localize, each one in a different well. 
The energy of the eigenfunction in the well at position id 
(with i = 0, ±1, ±2, . . .) is approximately given by 

Ei^^h[;; + Ve.tm, (2) 

where Uoj/2 = y/sEr is the energy of the lowest level of 
the lattice well. The first important observation is that, 
contrary to the uniform system case, the localized states 
are here exact eigenstates of the system for finite values 
of s. The amount of localization of the eigenfunctions 
depends on s and the external potential. The second 
crucial point is that, by increasing s, states with higher 
initial energy evolve into states localized farther away 
from the origin as compared to those with lower initial 
energy. For the register loading procedure, this implies 
that the constraint tp = E^ oi the homogeneous system 
is relaxed to ep < Er, since the missing particles result 
in empty holes localized at the wings, and not randomly 
distributed as before. These facts change dramatically 
the loading conditions and are the major results of the 
present work. 

For large s, the energy spectrum jSJ is characteristic 
of all states in the first pseudoband, while the states in 
the second pseudoband have energy approximately given 
by 'ihCj/2 + Vexti^d). One could have hoped that the 
condition Er > ep could be further relaxed, since the 
energy of the states in the second pseudoband progres- 
sively increases for increasing s, leaving more and more 
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localized states of the first pseudoband energetically al- 
lowed. This is, in fact, not the case as can be seen by 
looking at the behaviour of the levels of set (II) . These 
initially correspond to states in different pseudobands. 
Due to the localization effect, as soon as a gap opens, 
there is very little spatial overlap between eigenfunctions 
with similar energy belonging to different pseudobands. 
The small spatial overlap corresponds to a very small 
tunnelling between almost degenerate states. This fact 
is also clearly reflected in the structure of levels in Fig.^ 
where for increasing s the gaps at the avoided crossings 
very quickly tend to disappear. For our loading purposes 
this means that over the time scales of interest, the atoms 
will remain trapped in their pseudoband even if, by in- 
creasing s, states in lower pseudobands have much lower 
energy. Thus if initially ep > Er, the ground state of the 
final potential cannot be reached in realistic time scales 
and one will be left with more than one particle in the 
central wells. 

The level structure of the combined potential and the 
Pauli principle have crucial consequences also on the ef- 
fects of a finite temperature, T, and on the adiabaticity 
time condition. We here summarize all the properties in 
a detailed comparison with the uniform system case: 

1. For large but finite s the eigenstates of the Hamilto- 
nian can be well localized, depending on the particular 
shape of the external confinement. 

2. An inner part of the register virtually free of defects 
can be generated since: 

(i) The strict requirement ep — Er, fixing the initial 
number of atoms, is relaxed to ep < Er; 

(ii) In an adiabatic process at finite T, the occupation 
number of the final localized states is equal to that of 
the corresponding level of the confining potential. For 
the states at lower energy in current experiments this 
can differ from unity by as little as one part in 10^ 0- 
As long as T < Tp and f{E>Er)<t:l (with / being the 
Fermi thermal distribution), the states with occupancy 
significantly different form one are localized in a region 
~ ksT around the Fermi level. 

(iii) Strict adiabaticity is not required. Since the ex- 
citations induced by a potential varying in a time of the 
order of t are only created around the Fermi level in a 
shell of the order of h/r, one just has to enforce that the 
raising process be slow compared to the inverse of the 
Fermi energy. For instance in the case of a ID harmonic 
oscillator, the condition becomes r ^ h/ep = l/{Nuj), 
allowing a factor of N faster initialization than in the 
bosonic case. 

We remark that the points (ii) and (iii) are direct con- 
sequences of the fermionic nature of the atoms, and do 
not have a direct analogue in the Mott insulator phase 
of trapped bosons . 

To complete our analysis of the loading it is necessary 
to discuss more in detail the localization of the eigen- 
states. In Fig. 121 we plot the density profiles at T = 
for two values of s, with a number of atoms such that 
ep = Er- We see that the density becomes flat in the 
presence of a deep periodic potential [l^ . A uniform 



density distribution, however, does not imply the per- 
fect initialization of the quantum register. Comparing 
Fig. Eld) and Fig. |3fa) (corresponding to the same pa- 
rameters), we see that even if the density distribution is 
flat, the lowest eigenfunctions are not yet localized. 
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FIG. 3: Localization of the lowest eigenfunctions for s = 12 
for a harmonic confinement as in Fig. Q and for a V-shaped 
confinement V = \z — zo\Er/{2Tv'^d) (b). The V-shaped po- 
tential is chosen to have the same slope as the harmonic one 
at 2 = 4d. 



The degree of localization of the eigenfunctions at site 
£ arises from the interplay between the tunneling S and 
the energy difference between neighboring wells at po- 
sition id, given by AEg « d[dVext{z)/dz]z=M- In the 
large s limit, the eigenstates will be localized in the sin- 
gle wells, apart from a small contamination with flrst 
neighbors. The probability amplitude for the atom lo- 
calized at site £ to find itself in the site £ ± 1 is given 
by a ~ 5/{2AEg). This simple relation, whose validity 
was also checked numerically, implies that it is the first 
derivative of the potential to be important for localiza- 
tion. In the case of a harmonic potential then, the atoms 
are very well localized in the outer wells, while the cen- 
tral wells always suffer some tunnelling, even if small, 
due to the vanishing derivative. In a V-shaped potential 
Vextiz) — /iA/Vi3|z — zo| (with /zm the magnetic moment 
and Vi? a constant magnetic field gradient), all wells are 
instead equivalent in this respect, as shown in Fig. |2| 

Consequently, while on the one hand in a shallow con- 
fining potential a larger number of atoms can be loaded 
following the above scheme, since more levels are found 
below the recoil energy, on the other hand, once the lat- 
tice is raised, a steep potential is required to get a good 
localization of the eigenstates. A convenient procedure 
to obtain both a large number of particles and a good 
localization is, thus, to first raise the optical lattice up to 
its maximum depth in the shallow trap (determined e.g. 
by LOin or V-Bi„), and then to make the external confine- 
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ment steeper (determined by w/i„ or \/Bfi„) in order to 
ensure good localization of all populated levels. 

With a very steep final confinement, it might happen 
that states localized in sites away from the center of the 
trap become higher in energy than states in the second 
or higher pseudobands localized in the center of the trap. 
This happens if Vext{(-d) > <i^/sEr. This is not a prob- 
lem as long as states in different pseudobands are local- 
ized in different regions of space and the tunneling be- 
tween those states is negligible. The outermost atom in 
the register localized at £max — N/2 is degenerate with 
the state in the second pseudoband localized around site 
f , given by Vexti^'d) + hij = Vext(im,axd) ■ The con- 
dition for negligible tunneling is Imax — ^' ^ 1- The 
explicit expressions for harmonic and V-shaped confine- 
ment are respectively ^maa:-^' = (8/7r2)ysA^(wj„/w/i„)^ 
and l,nax - f = (8/3)ViiV(VB„/V%,„)2. 

We now consider, as specific examples, '^'^K atoms in a 
lattice with periodicity d = A/2 — 431.5 nm and ^Li 



atoms in a lattice with periodicity d = A/2 = 335.5 
nm. For both species, we take the magnetic moment 
equal to the Bohr magneton jiB- The first step of our 
procedure consists in raising the optical lattice to its fi- 
nal depth. The next step is to make the external con- 
finement tighter. The conditions fixing the strength 
of the final confinement are a good localization of the 
eigenstates (a < 10~^) and a negligible tunneling to 
higher bands [imax — S> 1). In Table I, we show 
the parameters required to trap a maximum number 
of atoms TV = Er/ep = 10'^. The final optical lattice 
depth is s = 30, corresponding to a tunneling parameter 
5 = Ej.- Lighter atoms require stronger confine- 
ments in order to be well localized. Moreover, as already 
discussed in Fig. El the V-shaped potential is more ef- 
ficient than the harmonic confinement, since the eigen- 
states are localized to the same degree throughout the 
register. 
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TABLE L Experimental parameters for the realization of a quantum register with ''"K and ''Li atoms, where we fix TV = 10'^ 
and s = 30. The parameter describes the localization around site I (which is site-dependent in the harmonic case). 



The extension of our proposal to 2D and 3D lattices 
is straightforward, as far as the relation ep < Er fix- 
ing the maximum achievable number of atoms, the adi- 
abaticity condition and the effect of finite temperature 
are concerned. However, the localization properties of 
the eigenstates have to be investigated carefully. A 
higher dimensional register can contain a larger number 
of atoms. For instance, in the case of cubic lattice and 
harmonic confinement, the problem scales in a simple 
way with dimensions. In 3D, defining Sho = ^/h/rriLU 



and Lu = (wxOjyUJz)^^^, we have N ~ {■K'^a1^/2d'^)^ , which 
can be easily of the order of 10^. 
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